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Abstract 

A procedure is developed for using Soft Collinear Effective Theory (SCET) to generate fully 
exclusive events, which can then be compared to data from collider experiments. We show that 
SCET smoothly interpolates between QCD for hard emissions, and the parton shower for soft 
emissions, while resumming all large logarithms. In SCET, logarithms are resummed using the 
renormalization group, instead of classical Sudakov factors, so subleading logarithms can be re- 
summed as well. In addition, all loop effects of QCD can be reproduced in SCET, which allows 
the effective theory to incorporate next-to-leading and higher-order effects. We also show through 
SCET that in the soft /collinear limit, successive branchings factorize, a fact which is essential to 
parton showers, and that the splitting functions of QCD are reproduced. Finally, combining these 
results, we present a example of an algorithm that incorporates the SCET results into an event 
generator which is systematically improvable. 
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1. INTRODUCTION 



To test a model of particle physics we must be able to describe the distribution of particles 
it predicts. Only a fully exclusive distribution is truly useful, since in order to obtain a 
meaningful comparison between theory and experiment, the theory must undergo a detector 
simulation and be subjected to the same cuts as the experimental data. However, high 
energy colliders produce events with thousands of particles in each event. There is no easy 
way to describe even the phase space for such complicated final states, and therefore Monte 
Carlo simulations have become essential for the analysis of every high energy experiment. 
Practically speaking, the current generation of Monte Carlo tools has been able to reproduce 
the standard model remarkably well. But with the onset of the LHC, new energy scales and 
new kinematical configurations, such as events with many high pt jets, may appear, and 
these tools will be pushed beyond their validity. So it is our task to help these simulations 
incorporate as much information from analytical calculations as possible, including loop 
corrections and the cancellation of infrared divergences when appropriate, while still having 
them produce exclusive events. 

The problem with current techniques is that a number of not necessarily good approxi- 
mations are forced by practical considerations. For example, it is unreasonable to calculate 
the analytic expression for a thousand-particle amplitude, and so event generators resort to 
the parton shower approximation. This approximation starts with an event resulting in, say, 
two quarks. These quarks then branch into quarks and gluons, and evolve down in energy 
until they hadronize at some low infrared (IR) scale. The branching is treated as a classical 
Markov process with emissions governed by splitting functions derived in the strict coUinear 
limit; thus any quantum mechanical interference effects are lost. A lot of work has gone into 
improving the results from Monte Carlo simulations. Mostly, it has been directed towards 
incorporating higher order QCD effects to improve the distributions in regions where the 
parton shower cannot be trusted. And, generally, quite good agreement with data has been 
achieved. However, it remains an extremely important open question to estimate the errors 
in these techniques, and to be able to reduce those errors systematically. It goes without 
saying that comparing theoretical predictions to data, or comparing the output of different 
Monte Carlo schemes, is not the ideal method of estimating error when searching for new 
physics. 

Part of the difficulty in simulating QCD is that it becomes strongly coupled at large 
distances. But even when QCD is weakly coupled, there is not an obvious perturbation 
expansion. Of course, we can expand a differential cross section as 

n 

But when there are multiple scales in the event, such as the relative momentum of pairs of 
final state partons, large logarithms may appear at any order. Typically, 

m<2n ^ 

where Pi is some kinematical variable, and Q is a fixed reference scale, such as the center- 
of-mass energy of the collision. Even if <^ 1 we may have as\og^{p/Q) ~ 1 and so 
the perturbation expansion breaks down. In this case, the large logarithms need to be 
resummed. Parton showers do this resummation, but without the right framework, it is 
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easy to lose track of which terms are accounted for and which are not. The right framework 
is an effective field theory which can sum the logarithms, incorporate finite corrections, 
and be extended with power corrections to fully reproduce QCD to any desired accuracy. 
Such an effective theory description was introduced in and in this work we will give more 
details and expand on that approach. 

An effective theory for event generation should not disturb the infrared properties of 
QCD, where strong dynamics and hadronization take place. Infrared divergences in QCD, 
such as a logarithmic dependence on the jet resolution scale, can be physical and must be 
reproduced. These divergences are either soft, for example, when a gluon's energy goes 
to zero, or coUinear, such as when a quark emits a possibly hard gluon in almost the same 
direction. The effective theory with all these properties is the Soft Collinear Effective Theory 
(SCET) 13, 1^, 1^, 1^- In this paper, we show how SCET can be implemented in an event 
generator so that the results are systematically improvable. 

We will demonstrate how at leading order, SCET is equivalent to the traditional parton 
shower approach. Parton showers write the probability for a branching as 



The splitting function represents the probability of a parton to split, and the Sudakov factor 
accounts for the fact that if a parton splits at a scale r, it should not have already split. We 
will see the splitting function derived from collinear emissions in SCET and the Sudakov 
factor reproduced from renormalization group (RG) evolution. One of the important facts of 
QCD, which gives rise to parton showers and which allows Monte Carlo techniques simulate 
multiple branchings, is that a given cross section factors into a product of probabilities, 
namely the probability to create an initial final state, multiplied by the probabilities to have 
subsequent splittings. This will rederived from SCET. 

Besides justifying the parton shower approach, SCET can be used to improve it. The 
problem of how to properly combine QCD matrix elements with parton showers is first cast 
in the language of scale separation. When the logarithms appearing in the differential cross 
section are large, the scales are widely separated, and the effective theory can be matched 
at the higher energy scale, and run down to the lower scale. The splitting functions used 
in the parton shower describe the long distance behavior of QCD (and of SCET), while the 
short distance physics is different in the two theories. Although the short distance physics 
of QCD is not fundamentally part of SCET, it can be fully reproduced through a consistent 
matching procedure. Thus SCET is valid at all scales, in contrast to QCD or the parton 
shower separately. 

We would like this paper to be self-contained and generally readable by both SCET 
and the Monte Carlo communities, as well as people familiar with neither field. So we 
attempt to incorporate a terse and incomplete review of event generators in Section |21 and 
of SCET at the beginning of Section HI and in Appendix |X| Section 0] gives a schematic 
presentation of the general idea of event generation in SCET. As much as possible, we 
sequester the detailed calculations to Section EJ which is the heart of the paper. The results 
of Section m are applied in various ways in Sectional We show first how Sudakov factors are 



da = A(r)P(r, z) , 



(3) 



where P(r, z) is a splitting function and A(r) a Sudakov factor 




(4) 
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reproduced from the renormalization group in SCET. A discussion is included of next-to- 
leading log resummation. We then discuss how the splitting functions and the factorization 
of successive branchings are understood from SCET. Next, we explore NLO effects, which 
include the cancellation of infrared divergences in physical observables. Finally, we use the 
SCET results to calculate some observables from parton-level results. We present the thrust 
distribution for 3-partons states, and the 2-jet fraction a2- Section IHl first summarizes how 
differential cross-sections are calculated, and then describes a simple algorithm that can use 
these cross section to distribute events in a Monte Carlo program. Finally we present our 
conclusions and outlook. We also include in Appendix ^ some kinematical relations and 
conventions that are used throughout the paper. 

2. INTRODUCTION TO PARTON SHOWERS AND EVENT GENERATORS 

In this section we review some aspects of how parton showers and event generators work, 
with an eye towards comparing to the SCET approach. This section contains no new infor- 
mation, and readers familiar with the subject can safely skip this section. 

An event is typically generated in three phases jE S H IM- First a simple hard process is 
selected at the parton level, with a probability proportional to its production cross section 
calculated using standard Feynman diagram methods. Second, the partons, which are taken 
to be highly off-shell at the hard scale, radiate additional partons and "evolve" down until 
their off-shellness reaches a low scale, which is typically chosen to be of order 1 GeV. Finally, 
the partons hadronize using some confinement model. We will be concerned here mainly 
with the second stage, the parton shower, and its interplay with the other two. 

Parton showers use classical evolution to describe the emissions of particles: they assign 
a probability for one particle to split in two. This probability depends on two kinematical 
variables. The first variable, r, is chosen to be a measure of the virtuality of the initial 
particle, and tends to zero if the two final particles are collinear. The second variable, 
measures the relative energy between the two final particles. Using that the splitting is 
independent of azimuthal angle, and that the branching probability is linearly divergent as 
r goes to zero, the probability is proportional to 



up to terms of order . Here, P(r, z) is our notation, while Pab is the traditional spin 
averaged splitting function, with a, h indexing the final state particles. For example, the 
quark splitting function in QCD is 



In order to calculate the classical probability of an initial particle with virtuality ri to 
branch at a particular value r2, one needs to include the probability that no branching 
has occurred at a larger value of r. This is similar to the well known case of nuclear beta 
decay, where this no-branching probability is responsible for the exponential decay rate. The 
probability of no-emission from ri to T2 is then given by an integral, known as a Sudakov 




(5) 




(6) 



where Cp = |- 
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factor 

A(r„r,) = exp (- £ ^ J i,^M^P^(,)^ (7) 

Here, the scale fi at which is evaluated can depend on both r and z. 

As an example, consider a quark- ant iquark pair with virtuality Q. The probability that 
a branching; occurs clt cL SCcllc r is given by the probability that neither quark branched at a 
scale greater than r, times the probability for either of the two to branch at that scale. Thus 
the differential cross section to have three final state partons is given by the cross section 
for two final state particles, multiplied by a Sudakov factor and the sum of the two splitting 
functions for quark and ant iquark emission 

da-s = da2 A(g, r)2 [P(r, z) + P'(r, z)] . (8) 

where P' is the equivalent of P for the antiquark emission. 

Parton showers rely on two crucial assumptions. First, one neglects the interference 
between the emissions off the various particles, and second, the intermediate states are 
taken onshell when deriving the splitting functions. For example, in the emission of a gluon 
off a quark- antiquark pair, the square of the full matrix element in QCD is replaced by a 
sum of two independent emission graphs 



ffi + 



(9) 



Both of these assumptions can be justified in the limit r — * 0, which implies that the virtu- 
ality of the branching particle is small compared to its energy. In that limit the branching 
particle becomes almost on-shell and the interference between the two QCD diagrams be- 
comes sub dominant. 

Note that the two diagrams on the right hand side are not gauge invariant and thus not 
well-defined. The parton shower circumvents this problem by summing only over physical 
polarizations, which restores gauge invariance, ipso facto. The splitting functions are well 
defined because the residue of the pole is gauge invariant and does not get a contribution 
from interference. But the r° and higher order pieces depend on conventions. For example, 
take the square of the relative transverse momentum between the two final state particles, 
p^, as the measure of the virtuality r. The branching probability is then proportional to 
P(r, z) with r = p"^. Now consider the different choice t = t, where t denotes the invariant 
mass of the two final state particles t = {pi + pjY /Q"^. These two variables are related by 

P^= "^"'7i',y''''' -M.-l) + 0(«') (10) 

Using this to change variables in P{t, z), we reproduce the same pole terms as in P{pt, z), 
however the higher order terms differ. Thus we cannot define the higher order terms in 
splitting functions in a consistent way. 

From the discussion so far we have learned that parton showers give simple expressions 
for the emission of partons. These can easily be turned into powerful computer algorithms 
using Monte Carlo techniques, which can be used to generate final states with many partons. 
One starts with a cross section for a process with a limited number of particles in the final 
state. In practice, event generators often start with processes with only two final state 



5 



partons. The virtuality of these two partons is chosen to be comparable to the hard scale of 
the interaction Q, and all additional partons in the final state are generated by the classical 
probabilities of particles to split into two particles with lower virtuality. This can be cast in 
a Markov Chain process, evolving the system from high to low virtuality, adding additional 
particles through the splitting functions. 

Different parton shower algorithms all use the physics described above, but they differ in 
which next-to-leading order effects they incorporate. First, they use different choices for the 
evolution variable r 3? Ill- Second, they use different choices for the scale /i at which as{f^) 
is evaluated in ((7j). Finally, for each emission in a parton shower, z and r are distributed over 
phase space assuming that the partons are on shell. But subsequent branchings require the 
partons to be off shell, so there is an ambiguity in how to assign the kinematics All of 
these effects only contribute at next-to-leading order, but they can give rise to considerable 
differences between the different programs in certain cases. 

Another aspect of QCD that has to be taken into account carefully arises for the emission 
of a soft parton The analog of the Chudakov effect from QED is that, when 

integrated over azimuth, soft emissions at large angle are sensitive only to the aggregate 
color charge of all the contributing partons. The result is that large angle soft emissions are 
suppressed. This is generally incorporated into parton showers by either using an evolution 
variable which corresponds to angle (as in Herwig or by explicitly vetoing emissions 
which are not angular ordered (as in some versions of Pythia jl5l|). 

The assumptions of parton showers restrict their validity to regions of phase space where 
each successive emission has to have virtuality much smaller than any previous one, Q S> 
Ti ^ r2 ^ . . . . Regions of phase space with r Q are only correctly described by the full 
QCD matrix elements. Thus, only an event generator which combines both the full QCD 
matrix elements with the parton showers will yield predictions which are correct to a given 
fixed order in perturbation theory, while also summing the leading logarithms correctly and 
allowing for the additional production of partons with small transverse momenta. Several 
techniques have been developed which allow such a combination of QCD matrix elements 
with parton showers, including original work which has been imp lemented in Pythia and 
Herwig EtL [isl and some more recent developments jiol. I20[l2l[ l23l|. One popular 



example is the CKKW procedure [19|], which provides a way to combine tree level QCD 
matrix elements with parton shower evolution. The idea is to compute the exact tree-level 
QCD matrix element for an event, including interference, but with a lower cutoff Tq on the 
virtuality between any two particles. One then reconstructs the dominant diagram using the 
kx algorithm and reweights the event by multiplying with appropriate Sudakov factors. 
Parton showers are then added to these matrix elements, vetoing all emissions with r > tq. 
This avoids double counting between the emissions contained in the QCD matrix elements 
and those generated by the parton showers. This algorithm has been tested against data, 
and overall the procedure seems to work very well |22l ]. 

A truly accurate next-to-leading order calculation would also incorporate loop diagrams. 
To obtain the matrix elements at NLO requires calculating the one loop correction to the 
lowest order matrix element and combining it with the matrix element which describes the 
radiation of one additional parton. The difficulty in going to NLO is that both the virtual 
contributions at one loop and the real emission have infrared divergences, which only cancel 
when both diagrams are combined to calculate an infrared safe observable , and this 

cancellation is difficult y to encode numerically. MC@NLO [2^ provides one solution (see 
also BHHUilllHIll). It uses the fact that the first splitting in the parton shower 
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reproduces the IR divergence of the real QCD emission, and can thus be used to devise a 
subtraction from both the real and virtual diagrams that render both these contributions 
finite. This yields two separate matrix elements, which are both finite and can be used as 
starting conditions for a traditional parton shower. The initial results appear promising js^. 
One problem is that MC@NLO works at the level of cross sections, not matrix elements, 
so it runs up against the possibility of having negative weights. It is also not clear how to 
generalize the procedure to higher orders. 



3. JET DISTRIBUTIONS FROM SCET: SCHEMATICS 

As we have seen in the previous section, the traditional Monte Carlo method uses splitting 
functions and Sudakov factors to generate a fully showered event. An event generator 
typically starts with an underlying hard process calculated using matrix elements of the 
full theory. The splitting functions then generate additional partons from these simple final 
states. The Sudakov factor, which is the probability of no-branching, is included, and it 
resums the large logarithms at leading order. 

The traditional Monte Carlo method uses splitting functions and Sudakov factors to 
generate a fully showered event. Both splitting functions and Sudakov factors can be derived 
in the limit of small transverse momentum, and parton showers are only correct in the 
limit where each successive branching has pr much smaller that any previous branching: 
Q ^ Pt^ ^ Pt^ ■ ■ ■ ■ Thus, only events with widely separated momentum scales can be 
described using parton shower techniques. 

The occurrence of these widely separated scales makes it natural to reformulate the 
problem in the language of effective field theory. The appropriate effective theory for this 
problem is the soft-collinear effective theory (SCET), which is designed to reproduce exactly 
long distance physics in the limit of collinear or soft radiation. SCET is a simplified version 
of QCD where only collinear and soft degrees of freedom are kept and all others have been 
integrated out. Any calculation using effective field theories requires three steps. First, one 
needs to match the effective theory onto the underlying theory. Matching determines coeffi- 
cients in the effective theory such that at some short distance scale it exactly reproduces the 
underlying theory. This ensures that the effective theory below that scale contains the same 
information as the underlying theory. The scale for this matching calculation is typically 
chosen to coincide with the hard scale in the problem, for example the center of mass energy 
of the collision. This ensures that the resulting short distance coefficients do not contain 
any large logarithms. Second, the effective theory is evolved to the lower scales arising in 
the process one wants to describe. This is achieved technically using renormalization group 
(RG) evolution, which sums logarithmic terms of the ratio of the low to the high scale. 
Finally, the matrix elements of the operators in the effective theory are calculated at the 
low scale. 

As we will show in this paper, these steps naturally correspond to the ingredients in 
traditional parton showers mentioned above. The matching calculation encodes the hard 
underlying process which we want to study. The solution to the RG evolution gives rise 
to evolution kernels, which are equivalent to appropriate combinations of Sudakov factors. 
And the resulting matrix elements in SCET have the property that in the collinear limit 
their squares simplify to squares of simpler operators, multiplied by the splitting functions 
of QCD. 

The first step is to match the full theory onto the effective theory at the hard scale Q. To 
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match, we introduce operators in SCET and choose their coefficients such that the matrix 
elements in QCD and SCET are the same. The matching can be done perturbatively in as 
to any order. Once we have matched at the hard scale, we no longer need QCD. To be more 
specific, consider the process e"'"e~ hadrons. In the standard model, this is mediated by 
a current 

J = qTq. (11) 

where F = in the case of an intermediate photon or F = gvl^ + QaI^I^ if the Z boson is 
included. The current has a non-vanishing matrix element in many final states, and we can 
compute {J\pi ■ ■ -pn) = {0\J'\pi ■ ■ -Pn) for any number of partons. In SCET, the process 
is mediated by operators which are constructed out of the fundamental objects of SCET: 
soft and coUinear gluons, collinear quarks, and Wilson lines. Wilson lines are required to 
ensure gauge invariance of SCET, and each collinear field needs to be multiplied with Wilson 
lines in the appropriate SU(3) representation. For example, collinear fermions and collinear 
gluons always come in the combinations 

Xn = Wn^n, A>:, = ^Wn[n-D,D^]Wl (12) 

In the SCET literature is commonly written as B-^. 

To reproduce the production of two collinear back-to-back partons one requires an oper- 
ator in SCET which contains two collinear fields, one for each of the directions: 

Ot^^ = XnTXn. (13) 

Each operator comes with a set of labels, corresponding to the directions of its collinear 
fields. We write 

Q{ri) ^ Q(ni,n2) Q(ni,n2,ni,) _ _ _ ^^^^ 

where Uj are the label momenta. The labels identify the degrees of freedom which cannot 
change - all the hard degrees of freedom which could have changed the collinear momentum 
in the full theory are integrated out of SCET. Note that there can be many operators with 
the same labels but different tensor structure (for example, the operators O-^ and O^^ we 
will define later on). For notational simplicity, we will often omit the labels (n^) when there 
is no ambiguity. 

The matching between QCD and SCET implies that we want to choose coefficients Cj 
for these operators such that the sum over all the operators reproduces the current of QCD 

J = G2O2 + C3O3 + C4O4 + ■ ■ ■ . (15) 

Since the Wilson coefficients Cj only encode short distance physics, they are independent of 
the choice of states used in the calculation of matrix elements. So, to satisfy this equation, we 
can take matrix elements in convenient states, and build up the Cj systematically. First, we 
take matrix elements with two quarks in the final state, which gives the matching condition 

{J\qq) = G^mqq) , (16) 

as no other higher order operators have matrix elements in a 2-quark state. This allows us 
to determine C2. Next, we take matrix elements with two quarks and one gluon 

{J\qq9) = C2{02\qqg) + Cs{Os\qqg) , (17) 
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and allows us to determine C3. The Wilson coefficients Cj with j > 3 are determined 
analogously. In order to correctly match QCD with up to m well separated partons requires 
operators Oj with n < m in SCET. 

The next step is running. The matching determines the Wilson coefficients at the match- 
ing scale, Cj{fi = Q) and the running will allow us to obtain them at a lower scales. The 
Wilson coefficients at these two scales are related by 

c„(/i) = c„(g)n„(g,/i). (is) 

The calculation of the evolution kernel n„ is a straightforward application of the renormal- 
ization group, and involves calculating the anomalous dimensions of the operators in SCET. 
Because the interactions in SCET are simpler than in QCD, the anomalous dimensions are 
fairly easy to compute. In fact, we will compute a closed form, algebraic expression, for the 
LL anomalous dimension of any Oj. 

The interactions in SCET allow for collinear gluons to be radiated off coUinear quarks 
and gluons, or collinear gluons to split into two collinear quarks. However, since SCET 
only describes the collinear or soft limit of QCD, the transverse momentum between the 
resulting particles has to be small. How small depends on the renormalization scale at 
which the emission is calculated, and the requirement is typically p± ^ /i. As an example 
consider the matrix element of O2 in an three parton final state with specific momenta 
\qqg)- Let px be the transverse momentum of the gluon with respect to a quark. ^ If the 
renormalization scale satisfies /i > pt, then gluon emission in SCET can give rise to a gluon 
in the final state, and a non-vanishing matrix element. If, however, fi < pt, then the SCET 
emission is not able to produce a gluon in the final state, and the matrix element of O2 
vanishes. 

For theory to be continuous across px, matrix elements at the scale fi = px — e must be 
equal to matrix elements at fi = px + ^- More precisely, we require that 

[c^m + c,{o,) + ■ ■ -i^+e = [C2m + c^m + ■ (19) 

To compensate for the discontinuity of (O2) at /i = px, we only need to change 03(03); 
all other matrix elements are continuous across the threshold. So we derive the matching 
condition 

[C2(02)]„+. = [C3(Of . (20) 

We call this matching threshold matching and it is done wholly within SCET. Note that 

('2') /f?Tm. 

the operator O3 in Eq. (j^UJ is different from the operator O3 arising in the hard matching 
ai fi = Q. To distinguish these operators we add a superscript, which labels the number of 
partons present at the hard matching scale. For simplicity, we omit the superscript for the 
operators matched at Q: Oj = Oj'K 

Let us apply these results and calculate the SCET expressions for a three jet final state, 
with the transverse momentum of the emitted gluon given by px- After matching QCD onto 
SCET at the hard scale fi = Q we find 

(SCET)q = C2(g)(02) + C3{Q){03) = (QCD) (21) 



Our notation is that pj_ is SCET notation for the perpendicular momentum label on a collinear field, 
while Pt refers to the transverse momentum between two four-vectors. 
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Using the RG evolution, we can obtain this at a lower scale /i 

(SCET)^ = C2(g)n2(g,/i)(02) + CsiQ)MQ,fx){Os) . (22) 

Finally, at the fi = px 'we have to perform the threshold matching to obtain 

(SCET)p, = C2(Q)n2(Q,pT)(Of ) + C,iQ)UsiQ,PT)m . (23) 

To calculate a differential cross section, we need to square this matrix element and sum 
over final state spins and polarizations. Using the Feynman rules in SCET we will show 
later that emissions in SCET factorize. In particular, 

\{02\qq)fP{pT,z) (24) 

where z = Eg/ [Eg + Eg). The function P{pt, z) is equivalent to a splitting function of QCD, 
up to power corrections, in the coUinear limit. 

We can now show that SCET agrees with QCD for large pt and the parton shower for 
small Pt- First, consider the limit pt ~ Q- Then, Iln{Q,PT) ~ 1 up to higher order 
corrections. So Eq. (j23p reduces to Eq. (j2H) and we reproduce QCD. Next, take the limit 
Pt ^Q- Then we are in the coUinear limit, and so the matrix element of O2 is very similar 
to that of QCD. Thus, {O^) ~ 0, since is the difference between QCD and SCET. So 
da ~ C2{Q)Il.\{Of^)'^ . Now, the kernel 112, is equivalent to the Sudakov factor, and as 
shown in Eq. (j24j) . the matrix element reproduces a splitting function, therefore the SCET 
differential cross section in this limit reduces to the parton shower result (jHl). 



4. REQUIRED CALCULATIONS IN SCET 

In this section we present details of the SCET calculations required to implement the 
scheme from the previous section. As we have seen in the previous section, there are three 
steps in the SCET calculations required to obtain jet distributions. First, we need to calcu- 
late the matching from QCD to SCET at some hard scale n ^ Q. Second, the renormaliza- 
tion scale of the operators is lowered using the renormalization group evolution. Third, a 
threshold matching is required when the renormalization scale gets lowered past the trans- 
verse momentum of one of the partons in the final state. Each of these three steps will 
be addressed in its own subsection. We will try to keep results as general as possible, but 
sometimes it will be necessary to choose a particular example. In this section we assume 
that the reader is familiar with the basic idea of SCET. For a quick review of SCET we refer 
the reader to Appendix [X] and the original literature 0, 0, 0, • This is the most technical 
section of the paper, and for readers not interested in the details, we will briefly summarize 
the results obtained in this section. This will make it possible to skip this section and still 
be able to follow the rest of the paper. 

The SCET results at 0{as) are as follows. The full QCD current is reproduced in the 
effective theory by operators On, such that 

J = C2O2 + C3O3 + of + ■ ■ ■ , (25) 
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The operators are given by 

O2 = Xn'^Xi 
= 9sXn, 



' 2 nq-Pt 



Xuq 
Xnq 1 



and for the Wilson coefficients required to NLO we find 



c.(g)^i-^ 

c-m = 1 

The Wilson coefficients C„ satisfy the RG equation 

d 



H — CnijJ^) = 7n(/i)C„(/i) 



The anomalous dimensions are 

■2 3 



72 = — 



as 

71 



73 = 

vr 

7n = 

vr 



Ci. ( lo, 

n 



,2 



11 
~6 



^C^ + ^C^)loK^ + 5. 



The Wilson coefficients evolve through the kernels n„: 

c„(/ii) = c„,(/i2)n„(/i2,/ii^ 

For a general anomalous dimension 



TT 



7n = f^nlog^ + Br 



we can solve n„ explicitly 
nn(Q,/u) = exp 
= 1 - 

where 



6 



(26) 
(27) 

(28) 



(29) 

(30) 
(31) 

(32) 

(33) 

(34) 
(35) 

(36) 

(37) 



asiQ) 



a,(/i)y A 



2 

-^-Bn lot 



MQ) 



2 



r„iog^^ + 25„iog ^2 



(38) 



«s(/i) 



as{Q) 
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i + £«.(Q)iog^ 



2r2, 



Only the piece proportional to is the leading log resummation, the B^, 
of next-to-leading logs. 



(39) 

resums a subset 
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A. Matching from QCD to SCET 



The first step is matching from QCD to SCET. As discussed in Sectional the matching 
condition is 

J = C2O2 + C3O3 + C4O4 + ■ ■ ■ , (40) 

and we want to choose Wilson coefficients in the effective theory such that QCD is reproduced 
at the scale Q to a given order in perturbation theory. QCD matrix elements with up to 
n well separated partons appear at 0{g^~'^) in perturbation theory. In order to correctly 
reproduce these matrix elements, we require operators in SCET with up to n coUinear fields. 
In the remainder of this section, we will explicitly perform the matching onto operators O2 
and Os, and comment on how to extend these calculations to operators with four or more 
collinear fields. 

As discussed above, each operator Oj depends on j labels for the directions of the j 
collinear fields. For each set of labels a different Wilson coefficient exists, so that in principle 
each product of operator and Wilson coefficient in Eq. (|4Up represents an infinite sum over 
the various labels. We will carefully treat the label dependence in the two-jet matching, but 
will then neglect the label dependence in the further discussions. 



1. Matching onto O2 at tree level 

The first step is to ensure that matrix elements with two quarks are correctly reproduced. 
Since the operators Oj with j > 2 have at least 3 collinear fields, they do not contribute to 

these matrix elements. Thus we need to choose Wilson coefficients C'^^ so that 

{J\qq) = J2^T''Hot'-'^\qq), (4i) 

rii 

where 

o?"'^^^ = Xn,rx„, (42) 

is a general basis of operators with 2 collinear fermion fields. To evaluate matrix elements 
on the right hand side of (|4ip. we need to know how jets Xn act on quark states The 
simplest prescription is 

Xn\q) = Sn^^n^^ , (43) 

where = p^/Eg. If J' produces two quarks, they must be back-to-back in the center- 
of-mass frame. So with this prescription, only the operators with n2 = ni get a non-zero 
coefficient. With these conventions, the matrix elements of the operator O2 is identical to 
the matrix element of the QCD current J' and we thus obtain 

C^'*^ = l, C^^'"^ = 0, rii^ns. (44) 

SCET fields also have labels corresponding to their pj_ momenta, but we choose to turn off 
any operator with p± ^ 0. 

Other choices for collinear fields acting on quark states are possible, and in fact (j43|) is 
by itself ambiguous. We can always write Xn as a jet in a different direction n' using ()A7|) 



Xn ^ 



1 + J^f_ 

n'-p 2 



Xn' , (45) 
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where p± is transverse to n'. Thus it would seem that Xn\(l) 7^ even if q is not ahgned 
with n. This ambiguity is resolved at higher order. Note that n-n' ~ p^, thus as long 
as n-n' -C 1 the two fields are equivalent up to power corrections. This implies that any 
operator XnJ^Xn2 with ni-n2 <^ 1 could be used as the operator Of course, in order to 
determine the matrix element of the operator O2 requires knowledge of how a SCET field 
with a given direction n creates or annihilates a final state with momentum in a direction 
n' . Any choice other than Eq. (|43|) will make the required calculations more difficult, so for 
simplicity we will stick to (jlHj) for the rest of this paper. 



2. Matching onto O3 at tree level 



Now that we have determined the matching onto the operator O2 we can proceed to 
match with 3-parton final states. The matching condition reads 



(Jlqqg) = C2{02\qqg) + CsiO^lqqg) . 



(46) 



The left hand side is given by the matrix elements for real emission in QCD from the quark 
and the antiquark leg 



(Jlqqg)' 



gsi'q 

-gs^/^' 



(Pq+PgY 



(47) 



For the right hand side, the operator O2 and its Wilson coefficient were determined in 
the previous section, but we need to evaluate the three-parton matrix element for 02- In 
principle, the additional gluon can be either coUinear or soft. We will first choose the gluon 
to be collinear and later check that the resulting matching condition is satisfied in the soft 
limit as well. The collinear emissions in SCET can come out of either a vertex from the 
Lagrangian or from a Wilson line. The matrix elements are extracted from the interaction 
vertices, and simphfy with the equations of motion i^^n = ^^t/i = 0. 



(02lqqg) 



g^ /„ 1 \2 >>n^f^ 



(02lqqg) 



(Pq+P 

n-{pq 



Pg) 



(Pq + PgY' 



£ r 



+ 



n-{Pq+Pg) n-pqn-{Pq+Pg 



n-pq 

n-Pq ' n-{pg+pg) n-{pq + pg)n-pq 



Tin 



(02lqqg) 



w„ 



g sin— Tin 



(02lqqg) 



Wo 



-gsin^—in 



(48) 



n-pg n-pg 

As the SCET Lagrangian is constructed to be gauge invariant, the matrix elements should 
satisfy a Ward identity. It is a straightforward check to see that (02lqqg)^ + (02lqqg)^ 
vanishes when A^^ ~ Pg- This holds separately for the quark and antiquark emissions, as 
it must since SCET is invariant under gauge transformations of the collinear fields in each 
direction separately. 

There are several issues that have to be resolved before we can use the above results to 
determine the matching of the operator O5. First, note that the collinear emission from 
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a Wilson line gives a divergence if the large label momentum {n-Pg for {02)^'' and n-pg 
for {02)^^) becomes small. In SCET, the coUinear phase space integration involves sum- 
ming over the large momentum labels, as well as integrating over the residual part of the 
momentum. However, in the sum over collinear momenta the value n ■ p = and p^ = 
has to be omitted [36[. Thus, the phase space integration will never reach a value for the 
label momentum such that the unphysical divergence is realized. However, keeping track 
of the condition n ■ p 7^ in the phase space integration can be difficult, and we therefore 
propose a scheme which will allow us to integrate the phase space naively, without having 
to worry about small label momenta. In this scheme one neglects the emissions from the 
collinear Wilson lines all together, and ensures gauge invariance by summing over only trans- 
verse polarizations when squaring the resulting matrix elements. By neglecting the {02)^ 
contributions, we subtract the divergences at small n-p. Any appropriate prescription will 
modify the effective theory only at short distances, and therefore differences in prescription 
can be absorbed by appropriately adjusting the Wilson coefficients in the matching condi- 
tion. In practice, however, we have to evaluate matrix elements at finite momenta, which 
are neither exactly at the hard scale where the matching is done, nor in the strict collinear 
limit, where the different prescriptions coincide. Therefore the physical results may differ. 
However, these differences are beyond leading order in the SCET expansion, i.e. they are 
power corrections. 

The second issue is that emissions in SCET cannot change the directions of the fermions. 
Since the operators O2 that were matched on in the two jet matching have the two fermions 
in back-to-back directions, the matrix element of O2 can only give rise to final state with 
back-to-back fermions. But if momentum is to be conserved when we include the gluon, the 
directions of the quark and anti-quark momenta can never satisfy Ug = exactly. This 
implies that in order to insist on PHj) . we need to change the direction of at least one of C,n 
or using (jl^ . otherwise we could never get a no n- vanishing matrix element with three 
partons in different directions. There are again various possibilities to deal with this, and all 
of them will lead to the same results up to power corrections. The choice we will adopt here 
is that directions of collinear fields which are not involved in the emission do not change. 
Then we rotate the emitting field into the direction of its final state. In other words, we use 

- L,^^^[r^n,l [LmC-n - [L,r]^^^en, , (49) 

where Q is some arbitrary operator. This equation shows how to replace a fermion in SCET 
with a fermion in QCD, which is necessary to compare matrix elements in the same external 
states. 

We now see that two different operators O^^'''^''^ and Q^^'''^''^ contribute to the same final 
state. Pictorially, 




The first will turns O2 into when the quark emits, and the second when the antiquark 
emits. We want to emphasize again that this is just a convention, and equivalent to any 
other convention up to power corrections. 

The transverse momenta are measured with respect to the directions n and n. This implies 
that when the quark emits, Pg +pf = and when the antiquark emits, pf +pf = 0. Using 
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these results together with ()49p . the quark emission and antiquark emission contributions 
to the matrix element of O2 can be simplified to 

io.m^^ ^ (50) 

Using the results of the SCET matrix elements given in Eq. ()50p and the QCD matrix 
elements given in Eq. ()47p. their difference reduces to the simple form^ 

{J\m9Y - mmg)'^' = 9s^j^ff^Lj^r^n, 

{Jlmr - Wqgg)'^ = -^s^^^Lr^A^n, (51) 

The final step in the matching is to choose a basis for the 3-jet operators Og"'''. The 
standard convention is to have C3 = 1 at tree level, so we take O3 = J ~ O^- Thus, we 
arrive at 

where was defined in Eq. (fT2j) . Note that O3 is the difference between the QCD and 
SCET emissions. Since the operator O2, together with the SCET emission of an additional 
gluon describes the infrared physics of QCD, the operator O3 only has contributions for 
large values of the transverse momentum. 

If the gluon is soft instead of coUinear, there are two emission diagrams in SCET. The 
matrix elements are 

Wqqg)'^ = 9sL^nn , Wqqgf^ = -gs'^^T^^^ ■ (53) 
n-pg n-pg 

This can be compared with the coUinear results given our convention of dealing with label 
momenta going to zero. In that scheme we only keep the matrix elements (O)^ and sum 
only over transverse polarizations of the gluons. Taking the soft limit of {02\qqg) (by taking 
Pg <^ n-pq,n-pq) we find 

{O^lqqg)''^ ^ gsL^^^n , {02\qqg)'"' ^ -9sL^^^n • (54) 
n-pg n-pg 

Thus, the soft limit of the coUinear Lagrangian emission is identical to the soft gluon emis- 
sion. However, care has to be taken about the second part of the convention, namely the fact 
that we only sum over transverse polarizations. Since for gauge invariant amplitudes the 
longitudinal polarizations do not contribute to physical processes, the soft gauge invariance 
of SCET is enough to ensure that the sum over transverse polarizations is equivalent to the 
sum over all polarizations. SCET is by construction gauge invariant, but one can also see 
the invariance trivially from Eq. ()54|1 . Replacing ^ Pg, the two contributions are equal 
and opposite in sign and thus cancel. Thus, the amplitude satisfies the Ward identity and 
is gauge invariant. With these two observations it is obvious that the coUinear emission 
reproduces the soft emission properly. 

In summary, we have chosen the conventions 
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• All fields appearing in operators have labels for collinear momenta only; in our basis, 
all p± momentum components are zero. 

• Emissions are calculated with the collinear emission from the SCET Lagrangian only. 
Wilson line emission and soft emission are discarded. 

• To maintain gauge invariance, we include only transverse polarizations of the gluons. 

• Fermions in SCET are rotated to a direction aligned with their 4-momentum before 
matrix elements are taken. I.e. ^„ — > ^■Cn,, where = -^q^ for massless fields. 

All other conventions are equivalent up to power corrections. 



3. Matching onto O4 at tree level 



We can outline the calculation for matching 4-jet operators O4,. These will make up for 
the difference between what SCET predicts when we have used only 3-parton QCD results, 
and the true 4-parton QCD prediction. Starting from gg, a 4-parton state can be either 
qggq or qqqq. We need to solve 



(j|4) = Y,Gt\oP\A) + cP{Of\A) + Cl"\0f |4) 

(n) 



(55) 



Let us take the \qqqq) state as an example. Then there are two QCD diagrams which 
contribute , 



(56) 



{J\qqqq) = ^ 
On the SCET side, the Co's can contribute 




{02\qqqq) 




+ 




as well as (9-? 



and of course Oy. 



{Oslqqqq) 



(Oilqqqq) 




(57) 



(58) 



(59) 



Using the results for O2 and O3 obtained previously, we can determine C4 and O4. In a 
similar way, the tree- level matching for n jets can be worked out. Note that since we have 
already set the conventions for evaluating matrix elements, there is no additional ambiguity 
when we match to operators with 4 or more collinear fields. 
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4- 2- jet matching at NLO 



So far, we have performed the matching calculations at tree level, and the normalization 
of the operators was chosen such that all Wilson coefficients are unity. In this subsection, 
we will determine the Wilson coefficient C2 at 0{as), which will be required to describe 
jet distributions at next-to-leading order (NLO). The matching condition is still given by 
Eq. (j4ip . but the matrix elements need now be evaluated at one loop. Loop diagrams are 
in general both IR and UV divergent, and we regulate IR divergences by adding quarks and 
gluon virtualities while using dimensional regularization with (i = 4 — 2e for the UV 
divergences. As always, the UV divergences are removed with counterterms, and the IR 
divergences cancel in the matching. For renormalization and we will use modified minimal 
subtraction, with fi as the MS renormalization scale. 

The one-loop QCD vertex correction [s^l is 



in 



; + 21og-log- 



2 log 



log 



(60) 



We draw the photon explicitly, because this graph depends on choosing F 
The SCET diagrams involve collinear gluons 0: 



Y 



An 



2 2 2 
-— - - + -log^ 



-log^ 



1^' 



+ 21og- 



vr 



(61) 



in 



and soft gluons: 



2 2 2 -Pi 
-— - - + - log— ^ 
e e /i^ 



log^ 



1^' 



2 log 



-Pq 



4+ — 



(62) 



An 



^log^ + log^^ 



n 
~2 
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(63) 



We have not assumed the quarks are back to back in computing any of these diagrams; 
Uq ■ Uq appears in the final result only through Q^. The SCET graphs do not depend on the 
choice of F. 

The sum of the three SCET diagrams gives 



A, 



(2) 
O 



asCp 



An 



A 
e 



, 2-Q' 

log 



2, -Q' 

-log — — 
e /i^ 

4 log- 



5n^ 
~6~ 



(64) 



Note how the ^logp^ pieces, which are UV and IR divergent, drop out of the sum. This is 
a consistency check, as such divergences cannot be removed with counterterms. 



17 



We can see explicitly that the entire dependence on the IR regulators and is the 

same in Aj and Aq , as expected. Taking the difference between the QCD and the SCET 
amplitude we find 



(2) 



Air 



3 
e 



-log T- 



log — 31og- 



TT 



(65) 



The 1/e UV divergences are canceled by counterterm contributions, and matching aX fi = Q, 
we get 



C2(Q) = 1 - 



An 



111 



37rz 



(66) 



To determine the counter terms of the operator O2 we need the wave-function renormal- 
izations as well. For the matching, these were not required, since they are exactly the same 
in QCD an in SCET 0. One finds 



At: 



log 



—p 



1 



(67) 



Combining this with the vertex diagram gives the result that the UV divergences in the 
full theory cancel, as is expected for a conserved current. We find that the renormalization 
constant in QCD is just 1, while in SCET 



An 



2 2, -u2 3 
- + -log-^ + - 



(68) 



5. 3- jet matching at NLO 



We can also match QCD onto 3-jet operators at NLO. We will not perform the calculation 
here, but merely outline what is required and make some qualitative comments. The 3-jet 
NLO matching involves calculating the difference between QCD diagrams such as 

Af = jr-l + • • • (69) 

and SCET diagrams, such as ^ 



A^a^ = — ' ® ^ — + ■ ■ • (70) 

There are also contributions from the same diagrams required for 73, below. As for the 2-jet 
operators at NLO, the QCD calculation and the SCET calculation are separately divergent, 
but the divergences all cancel in the difference. 

Note, however, that in this case, there is no reason to expect the structure of the operator 
which results to be the same as O3. So we must add a new operator to the theory whose 
Wilson coefficient starts at 0{as). This operator, when squared (or interfered with O3), will 
produce a differential distribution for 3-jet events which has a different shape than the tree 
level 3-parton distribution. In this way, all the shape changes from loop-corrections in QCD 
can be reproduced in SCET. We will return to this discussion in Section IFO 
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B. Running 



While the operators in SCET reproduce the long distance physics of QCD, their short 
distance behavior is completely different. In particular, the operator O2 is divergent in the 
UV, while the the full QCD current q'-f^q is not. The matching performed in the previous 
section ensured that at some particular scale the matrix elements in SCET reproduce the 
matrix elements of QCD exactly at a given order in perturbation theory. All the difference 
between the short distance behavior of QCD and SCET above the matching scale fi = Q is 
absorbed in the precise numerical value of the Wilson coefficients Cn{Q)- If the matching 
would have been performed at a different scale, the difference between QCD and SCET 
matrix elements would have been different, since the amount of short distance physics that 
has to be accounted for is different in that case. Thus, the Wilson coefficients C„ must 
depend on the value of the matching scale fi. Since fi is just a renormalization scale, it has 
no observable physical effect and our final answer should be independent of this scale. This 
implies that the matrix elements in the effective theory should satisfy a renormalization 
group (RG) equation 

[C2(/i)02(/i) + C3(/i)03(/i) + . . . + C„(/i)a(/i)] = . (71) 

Since each operator On contains a different number of labeled coUinear fields and interactions 
in the effective theory cannot change this number, each contribution has to separately satisfy 
the RG equation^ 

/i|-[C„(/i)a(/i)] = 0. (72) 

The n dependence enters the operator On only through its renormalization constant Zn- 
The anomalous dimension of an operator is defined as 

This allows us to write the RG equation in its final form 

Ai^C„(/i) = 7n(/i)C„(/i) . (74) 
This differential equation can be written as an integral equation 



n„(/i2,/ii) = — v = exp 



r -^m] ■ (75) 

J til 



We call Un the RG evolution kernel. It determines the change of a Wilson coefficient C„ as 
the scale changes. Thus, if we have calculated the Wilson coefficients C„ at one scale, we 
can use this kernel to obtain its value at any other scale. 



^ Beyond leading order, there may be mixing among operators. In this case the anomalous dimensions 
would be matrices, but we stick to the case without mixing for simplicity. 
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The anomalous dimensions will have the form 



7n = - 



r(i) 



TT 



TT 



r(2) 



1 ^^ 



TT 



TT 



5(2) + 



(76) 



The first term in brackets, is often called the cusp anomalous dimension. It multiplies an 
explicit, linear dependence on log/x. This term arises because of the double poles in the 
renormalization constants Z„, which in turn can be traced back to the fact that full QCD 
has overlapping soft and coUinear divergences. Such a term in the anomalous dimension is 
not problematic if no higher powers of logarithms appear, because it can be resumed. And 
in fact, it has been shown that at any order in perturbation theory the anomalous dimension 
contains at most a linear dependence on such a logarithm js^]- 

As a practical matter, it is helpful to have an explicit form for the evolution kernel (f73|l. 
At leading order in as 



7n(/i) 



-(r„ log — + 



TT 



(77) 



where r„ and i?„ do not depend on fi. The integral over /i can then be performed explicitly, 
and we find 



n„((5,/u) = exp 



57r 



/5o'«s(g)" 



log 



asiQ) 



+ 1- 



as{Q) 



2 ^ , 

— Bn log — ^ 



(78) 



From the integrated expression, it is easy to see that the i?„ piece is subleading to r„; the 
latter has an additional log enhancement. Moreover, when we change the reference scale 
then the coefficient of log 02/01 will shift, showing that there are additional contributions 
at NLL of the same order as 3^. Thus, only the cusp anomalous dimension is required for 
leading log resummation. 

We will now determine the anomalous dimensions 72 and 73 at 0{a2) and the cusp 
anomalous dimension for 7„. 



1. Calculating 72 

The counterterm of the operator O2 was already obtained in the previous section as a 
byproduct of the calculation of C2 at one loop. We found in Eq. (j68|) 



72 (/i) 



as(/i)Ci; 



TT 



1 -^'-.3 



(79) 



This has the form (fTUj) as expected. We can then plug into (fTHj) with 



12 = Cf, and B2 = Cf\- + tii] . 



(80) 



Note that the anomalous dimension is complex. This is due to a physical effect. Because 
degrees of freedom have been integrated out, some fields which would go on shell when we 
cut a QCD diagram are no longer around. But for the theories to be the same, the imaginary 
parts of the cuts must be made up for somewhere else, and so they show up in the anomalous 
dimension. 
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2. Calculating 73 



To get the anomalous dimension for O^-, we need to calculate its renormalization constant. 
We will again regulate the UV by dimensional regularization, and add a virtuality to all 
external particles to regulate the collinear and soft IR physics. The coUinear graphs involving 
the fermions are the same as for O^^ and given in (jblj) and The collinear graph involving 
the gluon is (sT 




ill 



- + - log - log — ^ 



log 



-Pa 



TT 

~6 



(81) 



The soft graph across the fermions is the same as (|63|) . up to group theory factors: 




2 2, 
+ - log- 



-PgP^ 



The soft graph between gluons and quarks are: 



log^ 



-pIpI 



lJ''^{Pq+Pqy 



71 

T 



(82) 




4:71 



4:71 



Paf 



1, 2 -pIp" 



Paf 



— - - log 



1 , 



TT 

' T 



(83) 



(84) 



e ' IJ^'^iPq+PgY 2 IJ^'^iPq+PgY 4 

Finally, we need the wave function renormalization of the collinear gluon, which is the same 
as in full QCD 



3 ^ 3 



-- + log- 



'PI 



(85) 



e 

where Uf is the number of flavors. 

Adding up the diagrams, the Mogp^ terms cancel, as they must. We find the renormal- 
ization constant for at one-loop is 



+ lo. 



1 + 7^ 

471 \ 6 e 

1 11 

^e^ Qe 

and the anomalous dimension is 

_ ,2 

Cf ( log 



3 2, -{Pg+Pgf 



as 

73 = 

71 



3 

{Pg+Pq?""^ 



^ ^^^ -iPq+P9)\Pq+P9f 
e fJ''^{Pq+Pq)^ 



9a Aog -/^'(Pg + Pg)' 

2 V {Pq+P9Y{Pq + '^ 



3e 
11 



6 



This is of the form (ffH|l . with 
1 



^3 



^^2" ^ — — 



^2^6 ^ g2(p,+P5)2 



6 



(86) 
(87) 

(88) 

(89) 



Again, for leading log resummation, only the cusp anomalous dimension r3 is relevant. 
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3. Calculating the leading contribution to 7^ 

As we have discussed earlier, only the cusp anomalous dimension is required for a LL 
resummation. The cusp anomalous dimension is the coefficient of the l/e^ counterterm and 
gets contributions from both soft and coUinear diagrams. However, for every l/e^ divergent 
term in a particular diagram there is a term proportional to (logp^)/^ which cannot be 
absorbed into a renormalization constant. As we found in the 2-jet matching at NLO, this 
term is canceled once the coUinear and soft contributions are added. For the cancellation to 
occur, the total l/e^ terms from coUinear diagrams have to be (-2) times the total term 
of the soft diagrams. This implies that the LL contribution to the anomalous dimension 7„ 
can be obtained from coUinear diagrams alone. 

CoUinear fields in different directions do not interact with each other at leading order in 
SCET, thus the one-loop diagrams from coUinear gluons only involve one coUinear direction 
at a time. For example, the coUinear diagrams required for the renormalization of O2 were a 
diagram involving only the coUinear quark (Eqs. and (|67p). and the antiquark (Eqs. (|62p 
and (|67|)). For the operator an additional diagrams involving the coUinear gluon were 
required (Eq. (jHTj) and (jHS))). Note the wavefunction graphs do not have poles. 

Combining these results, l/e^ poles in the counterterm of a general operator On with Uq 
quark fields and Ug gluon fields is 



1 + ?^ 
An 



UnCp + riqCA 



(90) 



which gives the anomalous dimension at LL order 



7n = - ^ [uqCp + UgCA] log ^ . (91) 
Therefore, the cusp anomalous dimension for any operator in SCET is simply 

r„ = + "^Ca . (92) 



C. Threshold matching 

Each coUinear field from which the operators On are constructed can produce additional 
particles through interactions described by the Lagrangian of SCET. This implies that an 
operator On can contribute to final states with more than n partons. For example, we have 
used in the matching from QCD onto the operator O^ that the 3-parton matrix element of 
the operator O2 is non-vanishing, and it gave the dominant contribution for small values of 
Pt- However, in SCET only emissions which keep all external and internal partons near their 
mass shell are included, with the amount of virtuality allowed depending on the resolution 
scale of the effective theory. Since virtuality and transverse momentum are related, we use 
that only emissions with < fi are allowed in SCET. This implies that if the scale fi gets 
lowered below the pt of one of the final partons in the final state, the original matrix element 
no longer contributes to the final state, and a threshold matching onto an operator with an 
additional coUinear field needs to be performed. The threshold matching condition for the 
operator O2 was already given in Eq. ((20)), and for a general operator oii^ 

P^\0ii%=,,^, = [cSi(0Si)].=PT-. (93) 
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We now work out explicitly the matching of O2 onto 0\ . 

The threshold matching of O2 onto O.^ occurs at a scale px- Actually, there are two 
relevant pr scales for the branching: the transverse momentum of the gluon with respect 
to the quark, p^, or to the antiquark, p^. However, since the SCET results are only valid 
in the limit px -C Q, the two pxs are interchangeable up to power corrections. So we use 
Pt = min(p^,p^). We might also consider p^, the transverse momentum the quark (or 
antiquark) with respect to the gluon. But if p^ is the smallest transverse momentum, we 
are well outside of the validity of SCET, and this kinematical configuration is taken care of 
by the original matching to QCD. In fact, taking pt = min(p^,p^,p|,) is equivalent up to 
power corrections, and this quantity, to which spherocity reduces for 3-parton kinematics, 
has the added property of being infrared safe. 

Above the scale fi = px, the 3-parton matrix element of O2 is given by the sum of the 
two terms in Eq. ()50|) 



nriVq + Pff) p _ ^ nq-{Pq+Pg) iq 



^n, (94) 



Note that the l/n-{p + q) terms in Eq. (|94p are due to the non- locality of the intermediate 
quark propagator. However, the propagator scales as Q/p^, and for ji < px the non-locality 
of the propagator is less than the resolution of the theory. Thus, we can think of the matching 
of O2 onto O3 at /i = as serving to keep the theory local by adding a label. 

(2) 

This matrix element can be reproduced by the matrix element of an operator O3 , defined 



as 



= 9sXnq 



Xn, , (95) 



where Xn is the same quark jet appearing in O2, and A^^ is a gluon jet, that is a collinear 
gluon field wrapped in Wilson lines, as defined in (jl2j) . To satisfy the threshold matching 
condition we then find that the Wilson coefficient C3 vanishes for jj, > px, and for ^ = px 
it is equal to the Wilson coefficient of the operator O2 

Cf\px) = C2{px), Cl{fi>px)=0. (96) 
Using previous results, we can write 

Cf (/i) = C2{Q)U2{Q,px)MPT,fiMPT - /i) (97) 

Threshold matching is just another way of saying that O2 emits a gluon a.t fi = px and 

(2) 

turns into C3 But it is important to understand the emission process as matching, so 
that it is improvable. For example, it would not be hard to do the threshold matching at 
next-to-leading order. This would involve calculating one-loop corrections to both the O2 
emissions and to O2 ■ The NLO matching would allow us to go beyond the strongly ordered 
limit Pt ^ ^ ■ ■ ■ ^ Pt> which the parton shower is restricted, to correctly describe 
configurations with p^ ^ ' ' ' ^ Px^Px^^ ^ ■ ■ ■ ^ Pt- 
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5. UNDERSTANDING THE SCET RESULTS 



In the previous section, we worked out in detail some matching coefficients and evolution 
kernels in SCET. The important results were summarized in the beginning of that section. 
In this section, we will show how those results can be combined to obtain differential jet 
distributions and show that we agree with traditional perturbative calculations. We concen- 
trate on the process e^e~ —>■ partons, whose kinematics are reviewed in Appendix^ If we 
work to order we can have at most three partons in the final state (quark, anti-quark and 
gluon), and we can thus obtain the differential decay rate da/ (dsdt). We begin by showing 
the Sudakov factors used in parton showers are reproduced in SCET, and discuss the NLL 
resummation. We then show that the splitting functions are reproduced in the collinear 
limit. As a corollary, we derive within SCET the classical factorization that parton showers 
assume. We then discuss NLO results, and compare to QCD. Finally, we display results 
from SCET for the thrust distribution and the 2-jet fraction. 



A. Sudakov factors from renormalization group evolution 



First, we compare the RG evolution kernels obtained in the effective theory to the Sudakov 
factors which arise in traditional parton showers. We will show that to leading log accuracy 
the evolution kernels coincide with traditional Sudakov factors. 

To start, suppose we just match from QCD to 2-jet operators. Then, according to ^^i . 
at n = Q, Cj'" = 1 and all other Wilson coefficients vanish. For an emission at a scale pr, 
we need to run C2 down to = pt- Using the 2-jet anomalous dimension we find 



C^"H/i) = C2(Q)n2(Q,/i) 



with the RG Kernel given by 



n2(Q,/i) = exp 



TT 



d/i 
fi' 



log 



-/i" 



(98) 



(99) 



This result can be compared with the expression for the Sudakov factor in traditional 
parton showers given in Eq. ((71). The precise form of the Sudakov factor depends on the 
choice of evolution variables used, and on the precise value of the argument of Og that is used. 
Most Sudakov factors use for the scale of an approximation to transverse momentum given 
by 

T = z{l - z)t ^ p'^T ■ (100) 

As an example, we will consider the so-called NLL Sudakov 0, 0, 2^, which also uses r 
as the evolution variable. The limits of phase space are y/r/Q < z < 1 — y/r/Q. Thus the 
Sudakov factor becomes 



exp 




;ioi) 



The z integral can be evaluated analytically. Substituting jj = Vt' 

■2 3 



Pt gives 



A^^^(g,^0 = exp 



(102) 
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FIG. 1: Errors on the RG kernel (Sudakov factor) 112 from next-to-leading log uncertainties. The 
light band comes from one NLL effect, varying the '^s term in the anomalous dimension, from 
< i?2 < 3. The dark band is from another NLL effect, the r^a^log term in 72, varying 
between —2 and 2. We normalize so that = 1 corresponds to 1 TeV. The LO RG kernel is the 
lowermost curve in the figure. 

Because the integral is dominated for small values of /i, we can drop the power law fi/fio 
terms, as is done in the literature. Comparing (jl(J2|) to we see that the 0{as) evolution 
kernel reproduces this Sudakov factor exactly. 

The 3/2 term in ()102|) gives rise to subleading logarithms after integrating over /i, hence 
the name NLL Sudakov factor. However, there are additional subleading terms which are 
not included consistently. For example, changing the reference scale hr at which the renor- 
malized is defined gives rise to subleading terms which are also NLL. The optimal value of 

cannot be determined to the order we are working. Moreover, different Sudakov factors, 
based on different evolution variables, give a different constant term In other words, 
the Sudakov factor only gives the leading logarithms reliably, and the 3/2 term may be 
dropped. 

To illustrate this point, we show in Figure ^ the 2-parton evolution kernel (Sudakov 
factor) 112 with various NLL effects included. The light band shows the effect of varying the 
3/2 term from to 3. The darker band shows the effect of adding a NLL factor proportional 
to F2a^/(7r^) log/i'^/Q^ to the integrand in Eq. (j^ . varying F2 between 2 and —2 Note 
that the effect of the 3/2 term and the NLL F2 term are comparable at small /i. Since small 
fi is precisely where the Sudakov factors become important, it will be important to include 
all the NLL resummation in the Sudakov factors consistently. Later on, we will explore the 
NLL effects on the thrust distribution for 3-parton events (see Figure |2I). 

For consistency we will therefore work only at leading- log accuracy, and drop the 3/2 
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term. At LL all acceptable definitions of the Sudakov factors will give 

A];\Q,f,) = exp|^^'^^a.(/iOlog^} , (103) 

which is the LL SCET prediction as well. Analogously, the gluon Sudakov factor (the 
probability for a gluon not to branch) is given by [1^ 

In SCET, the evolution kernel n„ for an operator with n coUinear fields is given by 

n„(g,/i) =exp|-^'^^7n(/^')} , (105) 

where 7n(/^) is the anomalous dimension for the operator, given in Eq. (jH^ . Comparing 
Eq. (ITnK|l and ^ with Eqs. (ITnH|l and (ITni|) . we find that at LL accuracy the RG evolution 
kernel and the Sudakov factors are related according to 

n2(g,/i)=A^'(g,/i)A^«(Q,/i), (loe) 

where Uq [ug) are the number of collinear quark (gluon) fields in the operator On- Thus 
SCET reproduces all the classical Sudakov no-branching probabilities at leading log through 
the renormalization group flow of an effective theory. 



B. Splitting functions from collinear emissions 



In this section we will show how the couplings in SCET, when put into cross sections, 
reproduce the splitting functions of QCD. 

For a 3-parton final state, once ii has run below Pt-, 2-jet operator can no-longer contribute 
and the threshold matching turns O2 into 0\ . If we are not concerned with additional 
emissions, the differential cross section for emission will be 



d(T 



dsdt 647r2 



phys pols 



where 



(107) 



(108) 



We have already seen that |C3^^p 



matrix elements. Using the explicit form of O^^"' given in Eq. 
explicitly. According to the conventions of Section HJ we find 



encodes the Sudakov factor, so now let us look at the 
(2) -p^ ii>)iii\ ^g^^ perform the sum 



(Of) =8g-tCF 



s u 

tT7 



t u' 

~s~{s + ty 



+ 



it + u){s + u){s + ty 



(109) 



The first term comes from the square of the diagram with the quark emitting, the second 
from the square of the diagram with the antiquark emitting, and the third from interference. 
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Rewriting the amplitude in terms of t and z, and taking the hmit where the gluon becomes 
coUinear with the quark, so — ^ and t 0, the amphtude approaches 



da (To 



ds dt QA-K' 



(Of)f = ..Hi^il±^ + ... (110) 



where the ■ ■ • are higher order in Pt/Q- So we reproduce the QCD sphtting function (jHI), 
as required. 

Note that the interference term (the third term in Eq. (jlU9p ) does not have an s or t 
pole, so it is finite as — and represents a pure power correction. Since the interference 
is higher order in the SCET expansion, we may simply drop it. Recall that dropping 
interference terms is one of the approximations used in the parton shower, and so we see 
that it is justified by SCET. However, the interference term should be included following 
our conventions for evaluating matrix elements - if we dropped it, we would not reproduce 
QCD at the hard scale. Because we match at the matrix element level, it is important to 
keep the interference terms in. Note also that at leading order, we do not need to distinguish 
pip from pip. We only know that p%^p% ^ Q, which is where the SCET amplitude can be 
trusted. 

We can also show that another element of parton showers, namely that successive branch- 
ings factorize and may be treated classically, can be justified within SCET. Consider a 
general operator 



= Xnfi= 1>9 (111) 

where f2 contains any additional coUinear fields and the Dirac structure of the operator. 
For example, for the operator O2 we would have = Vxn- Now consider the emission of a 
coUinear gluon off the coUinear fermion The amplitude for this process is given by 



where the emission S'^ can be obtained from the Feynman rules of SCET. Explicitly, 



A 



nq -Pg riq- {pg + Pg) \ " Uq ■ Pq 



(113) 



The first term comes from the Wilson line emission^ and the second from the vertex in 
the SCET Lagrangian. An important property of S'^ is that its square has trivial Dirac 
structure. In fact. 



SlS'^ = 2^,^C^ ^"^^' X id4.4 - 8g',CF-l\^ X id4.4 • (114) 

stu Z t 1 ~ " 



The arrow represents the coUinear t — > limit, where we see that QCD splitting function 
P{t,z) appears (cf. Eq. ©), with an extra factor of 1/z. Now, consider squaring the 



^ We include the Wilson line emission here, in contrast to our previous conventions, to insure gauge invari- 
ance. The same results hold if we ignore the Wilson line but only sum over physical polarizations. 
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amplitude A and summing over spins. The spin sum gives a factor of "^q, which in the 
coUinear hmit is = zt^, which commutes with S'^. So, 

|Ap = Tr {^S^'QQ^Sl} = zTt {SlS^'iQQ'^} 

Since Sjj^S'^ is proportional to the identity matrix in spinor space, we can pull it out of the 
trace to obtain the final result 

\A\^ = 6A7T^P{t,z)TT{00^} (115) 

This shows that, in the coUinear limit, the amplitude for a quark to branch into a quark 
and a gluon is independent of the other fields in the process, and that the probability for 
branching is given by a splitting function. 

Having considered a single emission of a gluon, we can go one step further and allow 
for multiple emissions. If we call pi^^ the transverse momentum of the emitted gluon with 
respect to its mother particle, the multiple emissions can be treated as a succession of 
threshold matchings if 

# >p5^^^ > ... . (116) 

This ordering of transverse momenta is called strongly ordered limit. In that case, the first 
emission is encoded in the threshold matching from O2 onto , the next emission in the 
matching from onto , and so on. Since in the above calculation we have assumed 

a general operator Xn^, the results can be applied recursively to the square of the final 
(2) 

operator On that is obtained after all the threshold matchings are performed. Thus, the 
square of the final operator in the strongly ordered limit can be written as the square of 
the original operator O2, multiplied by products of splitting functions. This is precisely the 
result that a parton shower algorithm would give for the same process, proven within SCET. 



C. Comparison with full QCD perturbative results 

In this section, we will show how the SCET results produce cross sections which agree with 
QCD at next-to-leading order. We begin by calculating the most inclusive of quantities, the 
total cross section for e^e~ — > partons, at NLO. This will incorporate the NLO matching 
of C2 and the 2-jet anomalous dimension. It shows that all the 0{as) information from 
QCD is in fact contained in the effective theory. Then we consider the differential decay 
rate da{e+e- qqg), also at NLO QHIillElHIillilHIiJ. First, we reproduce 



the 1/e and l/e"^ divergences in this rate from the counterterms in SCET. Then we show 
that at 0{as) all of the large logarithms, which appear in the limit pt <C Q, are resummed. 

We begin with the total cross-section at NLO. This total cross section receives contribu- 
tions from both 2- and 3-parton final states, and the 2-parton states have to be calculated 
to one-loop accuracy, while for the 3-parton final states only tree level is required [2^ |2^ . 
In the effective theory this means that to work consistently at order we need the one-loop 
matching for C2 and the tree level matching for C3. We also need the one- loop matrix ele- 
ment for O2 and tree level matrix element for O^. Both the 2- and 3-parton cross sections 
(T2 and (73 are infrared divergent, but these infrared divergences cancel in the sum of the two 
terms. 
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We begin by calculating the 3-parton cross section a^. In SCET it is obtained by squar- 
ing the amplitude C2{02\qqg) + C2,{Oz\<i<ig) ■ Because of the matching condition pUj) . this 
amplitude is exactly equal to the QCD amplitude. Thus, regulating the IR divergence in 
dimensional regularization, the 3-jet rate is 



SCET 



(To 



2 3 



2 /i^ 



19 

Y 



in 

6 



(;i7) 



where we have written the result in terms of the 0(0°) 2-parton cross section with dimen- 
sionally regulated phase space 



[A-K 



3(1 -e)r(2-e) 
(3 - 2e)r(2 - 2e) ' 



:il8) 



For the 2-parton cross section we require both the Wilson coefficient C2 and the matrix 
element (0219^) at one loop. Since the 3-parton rate above was calculated for an arbitrary 
renormalization scale /i we need the Wilson coefficient at that scale. Combining Eqs. 
with '^.2{Q■,^J') we find 



C2(/i) = C2(Q)n2(g,/i) = i- 



47r 



log^ 



3 log 



;ii9) 



We also need the matrix element of the operator O2 at one loop. In pure dimensional 
regularization, the one-loop contribution to the bare matrix element vanishes, since in the 
effective theory all large scales have been removed from the theory and all infrared scales are 
set to zero. Thus, the only contribution to the matrix element comes from the counterterm, 
given in Eq. (jHHj) 



3 

e 



2 -/i^ 



(120) 



Combining these results, the 2-parton cross section is 

1 



a 



SCET 



a',\C2{02\qq)\' 



(To 



C2^{02\qq) 



bare 



(To 



2ti 



3 

e 



6 



8 + 2- 



+ log^^ + 31og^ 



where we have used = Z-^{On)^^'^. 

The sum of the 2-parton and 3-parton contributions is finite, and we find for the total 
cross section to order a. 



C^tot 



SCET 



SCET 



Air 



:i2ii 



which is the standard QCD result, reproduced in SCET. 

There is an easy way to see why (J2 came out the same as in QCD. In QCD, the vir- 
tual contribution to the 2-parton cross section is UV finite. Thus, in dim reg, all the e 
dependence comes from IR divergences. But the IR divergences in QCD are the same as 
in SCET. In SCET, no divergences appear at all when using dim reg, because there are 
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no scales in the problem. Equivalently, in SCET the UV and IR divergences precisely can- 
cel. So the IR divergences are equal to the UV divergences which can be extracted from 
the counterterm. Thus, the counterterm in SCET has all the information about the full 
dimensionally-regulated QCD answer, up to finite terms. And these finite terms are pre- 
cisely what is calculated in the matching. 

Just as SCET reproduces the virtual contribution to (J2 from QCD, it is also capable of 
reproducing the NLO contribution to da^. In QCD, this computation involves all 1-loop 
contributions to e"'"e~ —>■ qqg. Again, there are IR divergences, which are canceled when a^, 
the integral of tree-level 4-parton emission, and the 2-loop contribution to (T2 are added. In 
this paper, we have not included all of the relevant computations to reproduce at NLO 
completely, however we have enough information to reproduce all of the infrared divergences, 
as well as the dominant large logarithms in the finite part. 

The next-to- leading order QCD result for the dimensionally regularized 3-parton differ- 
ential cross section can be found in [s^li^]. Consider first the divergent terms. They are 
all proportional to the tree-level cross section 



dcr^ 



QCDj 



a. 



Ann 



2\ £ 



27r V 



(122) 



3C^ + ^-2C^log^ 



Ca log 



St 



+ 



The SCET prediction for the 3-parton cross section comes from 

d^scET ^ PS-3X |C2n2(g,PT)n3(pT,/i)(of |gg^7) + C3n3(Q,/x)(03|gg^?)|' 



PS^sX 



C2Il2{Q,VT)MPT,^^)^{Of\qqg) 

^3 



(123) 
2 



bare 



C,Il,{Q,^^) — mqqg) 
^3 



bare 



where PS| refers to the dimensionally regulated 3-parton phases space. The l/e^ and 1/e 
poles in QCD come from soft and/or coUinear IR divergences in the matrix elements. In 
pure dim reg, the bare matrix elements vanish in SCET, since the UV and IR divergences 
cancel. Therefore all the poles show up as UV divergences represented by the counterterm. 
Thus, at NLO, the poles are extracted from Eq. ()123|) by using Z3 at 0{as) and setting 
112 = 113 = C2 = C3 = 1- At order C(q!s), the counterterm, from Eq. (jSHl), is 



47r 



2C p -\- C / 



(124) 



+ - 



3Cf 



2 



2Cp log 



-u 



C^log 



-St 



{2Cp + CA)log 



As expected, \Z-^ reproduces all of the divergences of the QCD expression. Note that the 

(2) 

matrix elements for O3 and ^3 reproduce QCD, by the matching conditions, so we get the 
same daf factor in both cases. 

The other part of the QCD expression we should be able to reproduce are the dominant 
large logarithms. For the logs to be large, we need pt <^ Q- In this limit, SCET is a good 
approximation to QCD, and {Os\qqg) ~ 0. Then the large logs should be resummed in the 
Wilson coefficients, through the 112 factor. To order a^. 



n2(Q,/i) 



V05 + (3 



2m 



1 



:i25) 
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For prp <^ n = 



I stu 



[see Appendix \^ and the kinematical structure of the SCET 



cross section is that of a sphtting function. Thus, 



dafCET 



ds dt 



647r2 
aoP{t, z 



C2(Q)n2(g,PT)(02|gg^?)r + 



2tx 



stu 



stu\ 



log 7^ + 3 log 7^ 



Q 



(126) 



To compare to the QCD result, we need to extract all the relevant terms from the NLO 
QCD expression [s^. Since we can only reproduce the logarithmically enhanced terms at 
this order, we will work in the kinematic limit t -C s,m and expand the full QCD result 
in powers of t. First, there is a contribution from the pieces with the kinematics of the 
tree-level cross section. This includes the finite piece from Eq. ()122j) . evaluated at /i = stu, 
and an additional finite piece from js^. All the terms from these expressions with logt are: 



do". 



QCDi 



ds dt 



(To- 



asCp s^ + t^ + 2mQ2 



27r 



st 



(127) 
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stu 



stu 



stu 



Cf log^ — - 3 log — + log — + (2C^ - Ca) log — 



u 



Next, there is a part of the QCD expression which is not proportional to da^, but has a 
splitting function as its collinear limit. Its large logs are 



da. 



QCDa 



ds dt 



agCp + (s + u) 

^TT st 



t V 

-(2C^-C^)log^log^ 



(128) 



Finally, there is a contribution from the running of a^. The SCET expression is evaluated 
with as(/i) while the QCD expression with as{Q). Changing the scale for the QCD expression 
gives an extra factor: 



dcr 



QCD3 



ds dt 



a,(/i)CFs2 + t2 + 2MQ2 



2n 



st 



as (3o , stu 

2^T^°^e 



(129) 



Taking the t — s> limit of these three expressions and adding them we find 



do". 



QCD 



ds dt 



(To- 



asCp 1 1 + (1 - s) 



27r t 



271 



log' 



stu 

0" 



(130) 



which matches the SCET expression in Eq. fjl26|) . Thus SCET resums all the large loga- 
rithms. 

There are terms in the full QCD expression with logs of u/Q^. In configurations with 
u <^ Q either the quark and antiquark are back to back, or one of them soft. These 
configurations do not come from soft or collinear emission from O2, and thus the large logs 
of u/Q^ are not resummed by 112. However, large logs of u are resummed in U^; and once 
SCET incorporates all the 1-loop matching, all of these terms, as well as all of the finite 
terms in the QCD expression, should be accounted for. 
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D. Interpolating between QCD and parton-showers 



A useful way to explore the SCET prediction is by looking at infrared safe observables. 
We will compare QCD, parton showers, and SCET. Since we only worked out results for e+e~ 
to 2 or 3 partons, we can only compute observables that depend on the 3-parton differential 
distribution. Once SCET is incorporated into a full event generator, more complicated events 
can be produced. But the results of this section are sufficient to show, at least conceptually, 
how SCET interpolates between QCD and the parton shower. So we take our observables 
to be functions of the 3-parton kinematics: f{s,t). For example, thrust is 

T(s,t) = l--^min{s,t,l-s-t} (131) 

Now, let us consider the 3 cases in turn. 

First, take tree-level QCD prediction. Thrust is computed as 

1 da^c^ a^Cp f ^ ^ s^ + e + 2uQ 



(To dT 2tt 



[ dsdt^^i^-t^^5(T - T{s, t)) (132) 
J st 



The prediction for QCD at NLO would be the same, but normalized to the NLO total cross 
section ao{l + ^) instead of just ctq. 

The parton shower prediction, which resums the leading logs, is given by 



PS 



1 da 

(To dT 



J dsdt [A(g, p^)'P(t, z) + A(Q, 4)2p(s, z')] 5iT - T{s, t)) (133) 



where P{t, z) is the splitting function from the quark emission, and P(s, z') is the splitting 
function from the antiquark emission (here, z = and z' = ^2^)- 
Finally, the SCET prediction is 

dsdt\C2U2{Q, Pt) {Of\qqg) + C3n3(g, pr) mqqg)\'5{T - T{s, t)) 

(134) 



(To dT 647r2 

To evaluate this we use the evolution kernels 112 and from Eq. ()38|) with Eq. ()33l34j) and 
the matrix elements. The matrix element for O^^ is given in Eq. ()109|) : 



2 



st{s + t)^ {t + u){s + u){s + t)^ _ 

The others (evaluated with the conventions of Section Ej) are: 



(135) 



KO3l™>r = 169fC.|^^^g!l±0^ (136) 

mO?)(Olm + hx, ^ -^e,lC,0^±f>g^ . (137) 

Note that the only term that has infrared s or t poles is Eq. ()135|) . The others are finite, 
and not genuine predictions of SCET. They depend only on our conventions. 
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FIG. 2: Thrust distribution from 3-parton states, at Eqm = 1 TeV. QCD (dashed red line), parton- 
shower approximation (dotted blue line), and SCET (solid black line) are shown. The grey band 
is a representation of NLL uncertainties, by varying the and terms in 72 and 73, between 
and their true values. 



If it were possible to choose consistent conventions so that {0^\qqg) = 0, then only 

{O^^lqqg) would contribute to the distribution, and because of the matching, it would be 
the same as QCD. In this case, our distribution would be a Sudakov factor, multiplying 
the QCD cross section, which is just the prescription of CKKW [3]. So our prediction, at 
this order, is equivalent to theirs, up to power corrections. Nevertheless, we cannot simply 
choose (Oslqqg) = 0, because we need a consistent set of conventions which allow us to go 
to higher orders. 

The thrust distribution for these three approaches is shown in Figure El We also include 
in this figure the NLL uncertainties as a grey band. This band corresponds to varying the 
B2 and i?3 terms in the anomalous dimensions between and their true values. Although 
the true value is known, this variation represents the NLL uncertainty from the F2 and F3 
terms, which multiply a^log/i in 72 and 73, and which is not known. Confer Figure Q as 
well. 

There are two important features of Figure |21 worth observing. First, note how SCET 
interpolates between QCD and PS. Small thrust, which corresponds to large pt, is populated 
by events with hard jets. In this region, we expect QCD to be a good approximation, as 
there are no large logarithms, and PS to be a bad approximation, as we are away from the 
coUinear limit where splitting functions are derived. Note that SCET approaches QCD in 
this region. In contrast, large thrust is determined by events with soft or coUinear partons. 
Here, tree-level QCD is inaccurate, because of large logarithms, while the PS is closer to 
reality. SCET matches the parton shower here. So SCET smoothly interpolates between 
QCD and PS. 

The second feature worth noting is the effect of the NLL resummation. Notice how the 
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FIG. 3: Percentage of events which have 2-jets, as a function of cutoff, using the kx algorithm. 
Shown is QCD (dashed red hne), the parton-shower approximation (dotted blue line), and SCET 
(solid black line). The grey band is the NLL uncertainty, as in Figure [21 



grey band is large for T ~ 1, which is where the PS is supposed to be accurate. Thus 
it seems the PS is vahd neither at small T, where hard emissions dominate, nor at large 
T, where NLL resummation is important. In contrast SCET is valid at small T, but also, 
once all NLL effects are consistently incorporated, it has the potential to be accurate in all 
regimes. This shows that resumming next-to-leading logs may be crucial to get accurate 
distributions. 

Next we compute the 2-jet rate. To do this we need a jet definition |4^. We will use the 
the fc-r (Durham) algorithm 2^1. It defines for any two partons a and b 



y''' = ^ mm{El Et}{l - cosOab) > Vent 



Q 



We then use 



y^mm{y''^,y^^,y'>^-} 



:i38) 



(139) 



For a 3-parton configuration, ii y > ycut, there are more than 2-jets, otherwise there are 
only 2. Thus we compute the 2-jet rate by integrating the distributions for QCD, SCET, 
and parton showers, over the appropriate range. We avoid the infrared singularities by 
integrating over events with kj- > y^ut- For example. 



a. 



QCD 



OtsCl. 



dsdt 



s'^ + + 2uQ' 

7t 



0(y - Vent) 



(140) 



The results for this observable are shown in Figure 01 The effect of the Sudakov suppres- 
sion can be seen on the left-side, at small ycnt- Here, the tree-level QCD prediction drops 
below zero, showing that its estimate of (T2 is no longer trustworthy. Both the PS and the 
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SCET curves are still positive at low energy, implying that they give a better estimate of 
the 2-jet cross section. For large ?/cut, the 2-jet rate should be accurately given by QCD. The 
PS gets the rate wrong, because it is integrating over hard emission where it is not valid. 
The finite difference between the PS and SCET curves comes from the integral over large 
Pt, where the PS cannot be trusted. 



6. TOWARDS SCET EVENT GENERATION 

We have shown in the previous section that SCET reproduces parton showers, Sudakov 
suppression, and NLO QCD results, and that it smoothly interpolates between the hard and 
soft regimes where QCD and parton showers are valid. In this section we will illustrate how 
the SCET formalism is naturally suited to implementation in an event generator. It has 
the capacity to produce particle distributions with several high pt jets, while summing the 
leading logarithms. And it is at least as powerful as parton showers for recursively adding 
additional soft or coUinear partons. 

What we are after is the differential cross section da for events with an arbitrary number 
of final states. In an event generator, this amounts to computing the weight, or probability, 
of an event, given the kinematics of the particles in the event. We will review schematically 
the analytical expressions required and give a simple algorithm to obtain events distributed 
according to these distributions. 



A. Analytical distributions 

Suppose we want the weight for an event with r particles from a production process 
at center of mass energy Q. To calculate this, the SCET approach requires a sequence of 
matching and running steps. First QCD is matched onto SCET at the scale fi = Q hy 
requiring that matrix elements with a given number of particles are correctly reproduced 
by the effective theory. In practice, how many particles we include depends on the number 
of matrix elements that it is feasible to compute in full QCD (and in SCET). This will in 
general be much less than r. Let m be the maximum number of partons which are matched. 
So the matching turns on operators O2 through at fi = Q. 

After the matching, each of the operators is evolved to lower scales using the renormal- 
ization group equations. To do this, we need the threshold matching scales, which can 
be derived from the event's momenta, for example with the algorithm 2J|. So we have 



Pt^ > p^'* > ■ ■ ■ > p^^\ At /i ~ p^"* a threshold matching is performed where the operator 

(2) 

O2 is matched onto an operator 0\ ■ The set of operators is then evolved to the scale 
jjL ~ at which scale the operators and O^"* are matched onto O^l'^ and 0^^\ This 
continues until the scale /x reaches the scale = Ajr. The differential cross section is then 
(up to phase space factors) 

m 

da, = l5^C(^)(/i)0(^)r (141) 

= |[C2(g)n2n3 ■ ■■iir\of^ + [C3(g)n3 ■ ■ -n^jof + . . . + [c„(Q)n^ ■ ■ -n^joff . 

Note that since only operators Oj with j < m were included in the matching onto SCET, 
this expression will only reproduce the full QCD results exactly with up to m partons in 



35 



the final state. All additional partons are described by radiation in SCET thus rely on an 
expansion in pt- As we have shown, after squaring the matrix elements, the SCET radiation 
is equivalent to the splitting functions used in parton showers. So, 

dcrr.>m = dcTm X Partou Shower (142) 

This distribution will be valid if the px satisfy 

p«, . . . » » . . . » pP (143) 

So if we want the distribution for r hard jets, we should aim for m > r. 

B. A sample algorithm 

Now let us consider how to incorporate these results into an event generator. This means 
we need an algorithm for sampling the phase space and unweighting events. We do not 
intend to present a complete or ideal solution, but rather sketch one possibility. Suppose 
that the we have QCD and SCET matrix elements for up to m partons, then an algorithm 
might look like 

1. Start with a hadronic event at center-of-mass energy Q, according to as accurate 
calculations of total cross sections cis cirG available. 

2. Pick a configuration of m partons with probability proportional to dam in Eq. 1)1411) . 

3. If the minimum px of the configuration, plp is more than than IR cutoff Air, then 
we have an event with m or more partons. 

• Start a parton shower at /x = pp~'^\ using the SCET RG kernels as Sudakov 
factors. 

4. Otherwise, we have an event with less than m partons. 

(a) Pick momenta for m —1 partons with probability proportional to dcTm-i in 

Eq. dnn). 

(b) If the minimum pt of this configuration, p^"^^ is more than the IR cutoff Air, 
then we have an event with m — 1 partons. 

• Do not shower. Go straight to hadronization with this kinematics. Or, 
equivalently, shower starting at /i = Air. 

(c) Otherwise, we have an event with less than m — 1 partons. Return to step (a) 
with m ^ m — 1. 

There are several points worth elaborating on here. First, note that the differential 
distributions in SCET do not diverge as two partons become coUinear or one becomes soft. 
This is because we are not using the QCD amplitude and correcting with Sudakov factors 
later on, but using the SCET amplitude which has the Sudakov suppression built in. Since 
the Sudakov factor vanishes exponentially in either of these limits, it overcomes the linear 
power divergence present in the full QCD amplitude. 
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Second, although we need to use some type of jet definition to sort the pt in the kine- 
matics, there will be only a very weak dependence on this definition in the final results. The 
strongest dependence on the jet definition should cancel from matching to the hadronization 
routine. The detailed sorting of pt only has an effect on momenta which are strongly or- 
dered, in which case all infrared-safe jet definitions should agree. There may be a subleading 
effect, which amounts to power corrections, and is beyond the order we are working. Also 
note that because the amplitude is finite even for small pt, we do not need to employ an 
intermediate jet cutoff scale in the theory, between Q and Air. We emphasize that this 
algorithm produces fully exclusive events, independent of jet definition. For observables 
with jets, the events can be combined using whatever algorithm is desired, and if the same 
algorithm is used on the data there should be good agreement. 

Third, if we are working to a consistent order in a^, the distribution in step 0] should 
really be done at NLO for m — 1 partons, and NNLO for m — 2 partons, etc.. As we have 
shown, all the NLO information from QCD can be reproduced in SCET, but one still needs 
to integrate over the regulated phase space for m partons in the singular region to cancel 
the NLO divergences for the m — 1 parton event. Quantitatively, the NLO effects are small 
for the vast majority of observables, so this algorithm as presented should be good for most 
practical purposes. 

Finally, the reason we do not just skip step HI is as follows. The distribution of m partons, 
from steps 1 through 3, is as accurate as possible with the information given. But because 
of the Sudakov suppression, the amplitude to produce fewer than m jets is probably very 
small. In fact, in steps 1 to 3, we make a Monte Carlo estimate of the m-or-more parton 
cross section. The distribution at small pt is also sensitive to subleading logarithms, and so 
we probably cannot trust it. Therefore, step E] recalculates the distribution from the hard 
scale, which should be a better estimate of the m — 1 parton event shape. 

7. CONCLUSION 

We have shown how to construct an event generator based on effective field theory. The 
correct effective theory reproducing all coUinear and soft divergences of QCD is SCET, and 
we have shown that SCET at leading order it is equivalent to a conventional parton shower. 
The advantage over the parton shower is that an effective theory is systematically improvable 
order-by-order. In particular, we have shown that by matching SCET onto QCD matrix 
elements with up to three partons allows to obtain cross sections which smoothly merge 
the parton shower approximation with more accurate matrix elements for events with large 
transverse momentum. 

Since we have a consistent effective field theory framework for the shower, we know pre- 
cisely to what order we are working. We can therefore estimate errors, and then incorporate 
higher order corrections if necessary when comparing to data. There are five places where 
higher order corrections can be incorporated 

1. Matching to QCD for higher multiplicity matrix elements. 

2. Matching to QCD at higher loop order. 

3. Higher loop running. 

4. Matching within SCET across the emission threshold, at higher order. 
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5. Power corrections from the SCET expansion. 

The first three we have aheady discussed at length in this paper. Of these the higher loop 
running may be the most important to work out and implement. We showed in Figures ^ 
and 121 an estimate of NLL resummation. It is clearly a large effect for small pt- Item 4. is 
not necessarily worth doing. We know that leading order SCET, like the parton shower, is 
best when p^"* 3> p^'* ^ ■ ■ ■ ^ . So if two successive emissions occur at similar pr's, the 
rate will be untrustworthy. However, a large pt emission will appear as another jet, and so 
we really should do the higher order matching to QCD, as in 1. and 2. to get the rate right. 
Item 5 is certainly important, but it may be prohibitively difficult to implement. 

Our results were collected in the last two sections, and we have indicated how SCET 
might be used to construct an event generator. We have also investigated the uncertainties 
in traditional parton showers. For large values of pt, these uncertainties arise from higher 
order in pr/Q, which are not included due to the approximations required to derive parton 
showers. For small values of pt on the other hand, the subleading logarithms which are not 
properly summed in the Sudakov factors give rise to large uncertainties as well. So, parton 
showers have relatively large uncertainties in most regions of phase space. Thus, while 
parton showers are very useful to populate phase space of events with many particles, it is 
very difficult to obtain results with small well defined theoretical uncertainties. Since SCET 
is able to improve the precision in all regions of phase space, we showed how SCET might be 
used to produce events with controllable errors. Our precscription was to produce an initial 
m-parton distribution at high accuracy, including resummation of large logarithms within 
SCET, and then to use the parton shower and hadronization routines, which are constrained 
by unitary, to fill out the jets with particles. One can go further, for example by using the 
properties of SCET, the accuracy of the shower routines could be improved, in principle, 
to any accuracy desired. Since effective field theories resum logs and are systematically 
improvable they have the potential to greatly improve our theoretical understanding of 
particle distributions. 
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APPENDIX A: INTRODUCTION TO SCET 

SCET is an effective theory containing only soft and collinear degrees of freedom, which 
can propagate over long distances. Since all long distance physics in massless QCD is 
determined by either soft or collinear particles, SCET is reproducing the long distance 
behavior of QCD, while all short distance physics can be encoded in short distance Wilson 
coefficients. SCET is essentially a simplified version of QCD in which all the IR degrees 
of freedom, that is, the soft and collinear fields, are the same, but the UV structure is 
simplified. In particular, fields which are neither soft nor collinear are integrated out. 
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The notation used in SCET emphasizes the nature of the SCET fields. Let be a 
hghthke direction = (1, rij), = 1. Any four- vector can be decomposed with respect to 
and = (1, — rij) as 

= ^{n-p)n'' + l-{n-p)n'' +p'^. (Al) 

So we can define 

ll^Y-^^n'^ + lin'', (A2) 
and derive useful identities, such as 

n-n = 2, 1^1^ = 0, {^,^}=4, {^,jf^^} = 0, (A3) 



p = {n-p){n-p) + pj_ . (A4) 

We also denote the normalized four-vector in the direction of p by Up. An operator V is 
often used to project out label momenta. For example, 

' -Xn=—^n. (A5) 



n-V n-p 

A field is coUiner to n if its momentum satisfies \p\_\ < Xn ■ p, where A is a small number 
giving the expansion parameter in SCET. The momentum of a coUinear field scales like 
{n- p,n- p,p±) ~ (A^,1,A), where we have used flA4j) with = 0. A field is soft^ if its 
momentum scales like p ~ (A^, A^, A^). Note that the sum of two collinear momenta in the 
same direction is collinear, so collinear interactions are allowed. But if two fields are collinear 
to different directions, rii and 77-2, they scale differently and interactions are forbidden. 

A Dirac fermion tpp with momentum p can be decomposed into collinear fermions ^„ and 

i^P = ^i^P + ^i'p = ^n + ^n. (A6) 

The difference betwen ipp and C,n is that on-shell, ipp satisfies "^p = while C,n satisfies 
ii'^n = 0. When n is aligned with p, ipp = but in general 7^ 0. The SCET Lagrangian 
is derived from the QCD Lagrangian by integrating out the small components H- Then, 
using the scaling properties of the quarks and soft and coUienar gluons, the Feynman rules 
are worked out as an expansion in A. They can be found in 0, 0,13,1^ • 

SCET has the curious property that no information is lost when C,n is integrated out, 
even if only the first order terms in the A expansion are kept. There are a few ways to see 
this. First, note that (^^)/4 and (^^)/4 are projectors, in that they are orthogonal and 
complete. However, we can regain ipp from just or using the identity 



^ Sometimes these fields are called ultrasoft, with soft denoting p ~ (A, A, A). These "soft" fields are not 
relevant for the current considerations as they cannot interact with either collinear or ultrasoft fields. 
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and that i/jp = ^np- One should keep in mind that SCET is a boosted version of QCD 
and any interactions between colhnear fiels in the same directions are thus given by full 
QCD. Whenever there are fields colhnear to different directions, however, their interactions 
will be different from in QCD. It is in this situation that SCET is useful. 

Since colhnear fields in different directions do not interact with one another, there should 
be a separate colhnear gauge invariance for each direction n. To ensure this gauge invariance, 
however, requires that each colhnear fermion is multiplied by a colhnear Wilson line, making 
the resulting field gauge invariant by itself 

Xn = Wn^n. Wn = exp j C/,^ I . (A8) 

Each colhnear fermion is therefore wrapped in colhnear Wilson lines, and it is these colhnear 
"jets" which are the basic building blocks of operators in SCET. For example, an operator 
with two jets would be 

o(-'-)=Xn.rx.,, (A9) 

where F is some tensor structure. Because ni and n2 are different labels, no colhnear fields 
can couple to both jets. Colhnear gluons can couple one jet to itself, and soft gluons can be 
exchanged between jets. 

Diagrams in SCET are usually computed in MS dimensional regularization. It is also 
helpful to regulate IR divergences by giving the quarks and gluons small virtualities. These 
virtualities must drop out of physical calculations, and it is a helpful check on the theory 
to show that they do. Sample calculations can be found in Physical quark masses can 
easily be accounted for in SCET [i^, but for simplicity we take all fields to be massless. 



APPENDIX B: 3-PARTON KINEMATICS 



For the three parton final state, label the momenta of the quark, antiquark and gluon p^, 
and p^, repsectively. There are two independent invariants we can construct. However, 

it is convient to go between a number of kinematic variables. Let Q be the center-of-mass 

(COM) energy of e^e~. Then, we have 

t = (Pg+Pqf =Q^il-Xq) 

U = {pq+Pq)^=Q^{l-Xg). (Bl) 

The invariants satisfy s + t + u = Q"^ and Xg + Xg + Xg = 2. The Xg are half the energy 
of the corresponding particle, in the COM frame. We also define the four light-like vectors 
pointing in the directions of the three particles as Ug, rig and Ug. So 

Pq = ^^q^q ' = ^^q^"^ ' Pg = J^O^g ' (^2) 

For each particle, we will also need the light-like vectors pointing in the opposite direction 

ni = (l,ni) fii = (l,-ni) (B3) 

The scalar products between different Ui can be obained from ()B1|) . and we also have 

nffij = ni-rij , nt-rij = 2 — nt-rij . (B4) 
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We can also derive ^ 

n^-{Pq + Pg) = nq-{Pq+Pg) = Q (B5) 

and its permutations. 

In the COM frame, the transverse momentum of the Pg or pq with respect to Pq, Pq or pq 
with respect to Pg, and Pq or pg with respect to Pq are given, respectively, by 



^'stu y/stu q Vstu 



XqC^ XgQ XqC^ 



We will also define pt as 

Pt = niin(p^, p^, pf,) . (B7) 

Note that for px ^ Q, 

Finally, it is helpful when discussing splitting functions, to use variables defined with respect 
to only two partons. If we are considering the quark and gluon, then we can use the invariant 
mass of the pair, and the energy fraction of the gluon 

t = Q\l-Xq), z = —^. (B9) 

9 ' 1 

Some useful relations between these variables and the others are 

^5 = 1-^, xq = z{l + ^), Xg = {l + ^){1-Z). (BIO) 
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